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In the field of Laurent series Q((T)), we construct a subring 5S wliich lias a natural 
Pj . D-module structure. Identities of Bernoulli numbers and polynomials are obtained 

^ ■ from the algebraic structures of OS. 
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1 Introduction 



Let e'^ = Yj^qT^/H be the formal power series with coefficients in the field 
Q of rational numbers. In the field of Laurent series Q((T)), the series B = 
B(T) = T/(e^ — 1) is contained in the formal power series ring Q[[T]]. The 



^ \ 2-th Bernoulli number B\ of order n is defined by 



oo r>Vn-) 
■pn \ ^ i rpi 

We write B\ simply as Bi. The z-th Bernoulli polynomial B\ [X) of order 
n is defined by 

oo rii'^) ( v\ 
-on XT _ y^ -P^ \^ ) rpi 

i=o '-• 
in the power series ring Q[X][[T]], whose coefficients are in the polynomial 

ringQ[X]. 
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Identities of Bernoulli numbers are abundant ever since early days of their 
discovery and still of contemporary interests. Among diverse approaches to 
these identities, we would hke to explore their algebraic structures. These 
identities, often of order one, are better understood through Bernoulli numbers 
and polynomials of higher orders. A framework for such a phenomenon is built 
up by the vector spaces defined below. 

Definition 1 We define 03 to be the Q-vector subspace o/Q((T)) generated 
by Laurent series T"^'B"'{bT)e°''^ , where m is an integer, n is a non-negative 
integer, a G Q and 7^ 6 G Q. We define Q3o to be the Q-vector subspace of 
<B generated 6?/T'"B". 

The power series B"(6T)e""^ G ® is the exponential generating function of 
UBi [a/b]^ which is the main ingredient of our identities. A multiplication on 
B"(6T)e"^ by a power of T simply means a shifting of the generating function. 
These observations motivate our definition. 

Clearly ®o is a subring of Q((T)). In Proposition [3l we will show that 03 is 
a subring as well. Note that B~^ = T'^^e^ — T^^ G OS. Hence, for a fixed /c, 
the exponential generating function of the sequence {S{n,k)}n>Q of Stirling 
numbers of the second kind is contained in 03. Indeed 

71=0 "'• '^• 

Besides the Q[T, e"^]-module structure, the field Q((T)) has a natural left mod- 
ule structure over the Weyl algebra D := Q{T,d/dT). It is clear that 03 is a 
Q[T, e"^]-submodule of Q((T)). In Proposition HJ we will show that 03 is also 
a D-submodule of Q((T)). The proof is essentially due to Lucas [TO', p. 62]. In 
terms of our terminology, he showed that 03o is a D-submodule of Q((T)). 

Starting with elements of 03 and performing algebraic operations on them, 
we will provide methods to represent the result as a D-linear combination of 
elements of the form r™B(6T)e"^. Equating coefficients in the Laurent series 
of the representation, we obtain identities of Bernoulli numbers and Bernoulli 
polynomials. For instance, the relation 

B^ = ("l - 
gives rise to Euler's identity 




m— 1 



2m^ 



E ( 2, jB2iB2rn^2^ = -(2m + 1)52^ (1) 



1=1 



by equating the coefficients of T^™ for m > 1. See Relation [TO] for details. 
Such a viewpoint applies to many identities found in the literature including 



generalizations of Euler's identity [3|14fl5fl6] , an identity from Eisenstein se- 
ries [13], identities from special values of Zeta functions [5|6] . and others [1]. 
One may implement algebraic operations of *B into computer programs, which 
provide new identities elegant or involved. 

It is plausible that our algebraic view may be extended to other classical 
numbers. We leave the possibilities to the reader. 



2 Algebraic Structures 



Note that elements T™B" form a basis for the Q- vector space *Bo. This fact is 
a consequence of the transcendence of B over Q[T], which is equivalent to the 
transcendence of e^ . In the Q[T, e"^]-module 53, elements T™'B"(6T)e""^ may 
have non-trivial relations. For instance, 

Be^ = T + B. 

See Section [3] for more examples. 

The ring structure of 23 is based on the following version of partial fraction 
decompositions. 

Lemma 2 For any positive integers ui,- ■ ■ , ra^, fci, ■ ■ ■ , kg, there exist Qi G 
Q[X] and positive integers mi, ii with ii < ni + ■ ■ ■ + Hs such that 

1 ^ ft 



(x«=i - 1)^1 • • • (x'^^ - 1)"^ Y (^™' - ^y 



PROOF. We may assume that ki,--- ,ks are distinct. The case s = 1 is 
trivial. Assume that s = 2. Let k be the greatest common divisor of ki and /c2- 
Since A - 1, X^^^/^^-^ + ■ ■ ■ + A + 1 and A^^^/fc)-! + . . . + x + 1 are relatively 
prime, there exist gQ,gl,g2 G Q[A] such that 



1 



9*0 ^ 9i ^ 52 



Let go = g*oiX'^), gi = (X^ - 1)"^^K^'), 92 = {X'^ - 1)"^^2*(^')- Then 
1 9o , 9i 92 



(X^'^ — l)"l(A^2 _ Uri'z fX'^ — l)"l+"2 (x'^'i — l)"i (X^^ — 1)"2 ' 

The lemma follows from induction on s. 



Since the denominators of ([2]) has a common factor X — 1, the decomposition 
in Lemma [2] is not unique. For instance, 



1 1 - iX iX 



1 _i 

2"" I 2"" _ 2 I 2 



(X-1)(X2-1) (X-l)2 X2-1 (X-l)2 ' X 

Proposition 3 ^ is a subring o/Q((T)). 



PROOF. We need to show that B''^{{ki/k)T)B''^{{k2/k)T) E 53, where 
ki, k2 are positive integers ad k is an integer. Choose polynomials Qi and pos- 
itive integers rrii, ii as in Lemma [2] for s = 2. Replacing X by e^^'' in ([2]) and 
multiplying both sides of the equation by {kiT)"'^{k2T)"'^/k"'^~^'^^, we obtain 
the required formula 



• A; '^ ^ A; '^ '^^ k ' ^ k ' ^rrii ^ k 

1. * 



The Weyl algebra Z^ = Q(T, d/dT) is generated by the multiplication operator 
T and the derivation d/dT as a Q-subalgebra of HomQ(Q((T)), Q((T))) . For 
f,gED, we write their product a.s f ■ g. The field of Laurent series Q((T)) is 
endowed with the natural left D-module structure. For instance, the module 
structure gives {d/dT)T = 1 G Q((T)). Do not be confused with the product 
(d/dT) ■ T, which satisfies the Leibniz rule 

d ^ ^ d 

r = 1 + T . 

dT dT 

The subalgebra of D generated by T is also denoted by Q[T], which makes no 
confusion with the subring of Q((T)) generated by T. For f,gE Q[T] C D, 
their product is also written as fg. This notation agrees with that of the 
module structure on Q((T)). Moreover, using the Leibniz rule, any element of 
D can be written as E /i ■ {d/dT)\ where fi e Q[T]. We write E fi ■ {d/dT)' 
also as J2fi{d/dTy without any confusion with the module structure. More 
notations: d'/dT^ stands for {d/dT)\ For F G Q((T)), we use t/F"/rfT" for 
{d/dTY¥. In our convention, dF^/dT^ = F. 

Proposition 4 ^o o,nd !B are D-submodules ofQ{{T)). As D -modules, *B is 
generated by r™B(6T)e""^ and Q3o is generated by T'^B. 



PROOF. For the first assertion, it suffices to show (rf/rfT)(B"(6T)e"^) G ® 
and ((i/(iT)B" G ^Sq. This follows from the formula 

J, (B"(6T)e'^^) = ((a + |)B-(6T) - n6B-(6T) - ^B"+i(6T)) e'^^, (3) 



which can be found in [TOl p. 62] for 6 = 1 and a = 0. The general case of ([3]) 
reduces easily to the case a = 0, which is obtained by dividing the equality 



nB"(6T)(e^^-l)" 
= n{bT)'' 

mr 

dT 
= rA(B"(6T)(e^^-l) 

= n6TB"(6r)(e^^ - 1)" + nB"+V6r)(e^^ - 1)" + Tie^^ - 1)"— B"(6r) 

dT 

by T{e^^ — 1)", c.f. [8:, Prop. 1]. The second assertion follows from 

B-+\bT)e''^ =(l-bT + — - -4^] (B^(bT)e''^) , (4) 

\ n n dT j ^ ' 



which is equivalent to (j3]). 



In Section [3], we will provide two special cases of partial fraction decomposi- 
tions which suffice to calculate the product of elements in 03. 

The first case is a decomposition of 1/(X" — 1)(X'" — 1), where m ^ n. Let £ 
be the greatest common divisor of m and n. 

Definition 5 For i = 1, we define gm,n to be the polynomial with degree less 
than m — 1 such that 



(X - 1)(1 + X + ■ • ■ + X™-i)(l +X + --- + X"-i) 

-L gn,m gm,n 



mn(X-l) 1 + X + --- + X"-! l + X + .-. + X'"-!' 

In general, we define gn,m = gn/i,m/i{X'^) . 
For example, 



5'2,3 = -2' fi'3,2 = 3(^- 1)5 

92,, = -\. ^5,2 = \{2r^ -X^ + X-2)- 
93,5 = i(X - 1), (?5,3 = ^(-2^' + X^ - X - 3). 



We have a decomposition 
1 



+ 



9n,r, 



(X"-1)(X'»-1) mn(X^-l)2 X"-l X™-1' 

Note that (71^^ = 0. For n > 1 and £ = 1, we may describe gn,m = ^o + cii^ + 
■ ■ ■ + a„_2-^"^^ using the complex number v = e^^^^^^^ by 



^a,\ / 



Oi 



1 V 

1 t;2 



,71-2 



,2n-4 



\ 



yan-2) yi t'"^^ ■ ■ ■ i;("-i)("-2)y \^(^;("-i)™ _ i)-y 






The second case is a decomposition of 1/{X^ — 1)^{X"' — 1), where f is a divisor 
of n. 

Definition 6 We define /i^ „ and /{ „ to be the polynomials with degrees less 
than k and n — 1 respectively such that 

1 



uik) 



(k) 



+ 



Jl,n 



(X-l)'=(l + X + --- + X"-i) (X-1)^ 1 + X + --- + X"-!' 

We define hf^ = h[%,{X^) and /fj = /|5/,(X0. 
For example, 



5 



c{2) 



1 



f^'i = -(2X3 ^ 3j^2 ^ 3j^ ^ 2). 
5 

We have a decomposition 



/i' 



(fc) 



(fc) 






(X^ _ l)fc(X'* _ 1) (X^-l)fe+i ' X"-l 



Note that h[l = oq + ai{X — 1) + ■ ■ ■ + afc_i(X — 1)*^ ^ can be described 
inductively by gq = 1/n and 

for i > 1 with our convention ("j =0 for n < i. Note also that 

Jl,n /j^ _ ^\k 



3 Relations and Identities 



We list some relations in 23 and demonstrate their implications on identities 
of Bernoulli numbers and polynomials. 

Relation 7 

B(-T)-r = B. 

It is straightforward to verify the above relation, which implies that Bi = —1/2 
and all other odd Bernoulli numbers of order one are zero. 

Relation 8 

e^B = B(-T). 

It is straightforward to verify the above relation, which implies the recurrence 
formula 



J:[^B^ = i-lrB^ 



1=0 \ 



for n > 0. Together with the vanishing property of odd Bernoulli numbers of 
order one, the above recurrence formula can be also written as 



" ^n 



E B. = B„ 



^=0 \* 



for n > 2. See [121 p. 18, Equation (4)]. 
Relation 9 

n-l 

J2 B{nT)e'^ = nB. 



i=0 



Relation |9] follows from the direct computation 

n-l rp nT _ 1 

Y, B{nT)e'^ = n-— — ^B(nT) = nB. 

U) e^ -1 riT 

It is equivalent to Y17=o Be'-'^^"^"^ = nB(-)e°-^, whose coefficients of T™ give 
rise to multiplication theorem 

n-l • 



p. 21, Equation (18)]. 
Relation 10 



B"+i(6T)e"^ =(l-bT + — - -^] (B"(6T)e''^) . 



This is equation (^. For 6=1, equating its coefficients of T\ we get the 
formula for lowering orders of Bernoulli polynomials 

Bt^'\a) = (1 - i)5f V) + (a - ^)-i?a(a) 

[121 p. 145, Equation (81)]. 

For a = and 6 = n = 1, we get Euler's identity 



i 



Y^ [ . \BiBm-i = (1 - m)Bm - m5„_i 

j=0 



from the coefficients of T™. Using the vanishing property for odd Bernoulli 
numbers of order one, we recast Euler's identity in the form of ([T]). 

For a = 0, 6 = 1 and n > 1, repeatedly using Relation [TOl we get gener- 
alizations of Euler's identity. See R. Sitaramachandrarao and B. Davis [T5] . 
A. Sankarayanan [H], W.-P. Zhang [16], and K. Dilcher [3j. 

The polynomial analogue [3 Proposition 2] 

J2 r\Bi{a)B^_i{h) = (1 - n)B^{a + h) + n{a + h- l)5„_i(a + h) 
of ([1]) also follows from direct computation: 



(Be'^^)(Be''^) = ({1 - T - T-^)b\ e('^+')^ 

= (l + (a + & - l)r - T A^ (Be("+^)^) . 

Relation 11 

-^ = T--UT,B), 
where 

n+l 

fn := (-1)" E(J - 1)! {Sin + l,j)U''-^^'-nS{n,j)U'^-^) V^ 

wt/i Stirling numbers S{n,j) of the second kind. 

Relation [TT] is given in |H1 Theorem 1] with /„ G Z[U, V] constructed induc- 
tively hj fo := V and 

for n > 0. Using the initial values S{n, 1) = S{n, n) = 1, S{n, n + 1) = and 
the recurrence relation 

S{n + l,j) = S{n,j - 1) +jS{nJ), 

one shows that the two definitions for /„ agree, of. [U Lemma 3.1]. 

Given m and n, 

for some ip & D hj Relation [101 In terms of Bernoulli numbers, this gives 
rise to identities described by T. Agoh and K. Dilcher [1]. For example, the 
identity 

"■ fn\ 1 1 

V . ]Bi+iBi+n-i = -{n- l)Bn - Bn+i - -{n + 3)Bn+2 
f^oVJ 6 6 

comes from the relation 

dBdB_/ 1 t/3 i£ T ^ ^ i\ 
'dT'dT~\"6 dT^~2dT^^6 dT ~ dT ~ 6 J ' 

Combining with Relation [TOl we have another relation 



9 



where 

^ \ 6 dT^ 2 dT^ 6 dT dT 6 J 

Equating coefficients of T^" with n > 4, we get the identity 

y.' (2n - 2)! B2i B2n-2i ^ {2n + l){n - d.) ^ 

^ (2i - 2)!(2n - 2i - 2)! 2i 2n - 2i 6n ^"' 



which was proved by H. Rademacher [13] using Eisenstein series and by M. Eie 
[S] using Zeta functions. See also [S]. 

Relation 12 

B(mr)B(nT) = B\iT) + mTB{nT)g^^^{e^) + nTB(mT)(?^,„(e^), 

where i is the greatest common divisor of distinct positive integers m, n and 
polynomials gm,n, gn,m o'^e given in Definition\^ 

For example, 



B(2T)B(3T) = B^ + -T{e^ - 1)B{3T) - -TB(2T), 

B(2T)B(5T) = B^ + -T(2e^^ - e^^ + e^ - 2)B(5T) - -rB(2T), 

5 2 

B(3T)B(5T) = B^ + -T(-2e^^ + e^^ - e^ - 3)B(5T) + -T(e^ - l)B(3r). 

5 3 

In terms of Bernoulli polynomials, the above expression for B(2T)B(3T) to- 
gether with the relation B^ = (1 — T — T{d/dT))B gives rise to 

1 3^2"- (^) B.5_ = |^5.-i(^) - (1^ + |i + -) Bn-. + (1 - n)5„. 
For even Bernoulli numbers, the identity reduces to 

^32.22n-2. j^2n^5^^^^^_^^ = ^i?2.-i(^) + (1 - 2n)i?2., 

where n > 1. 
Relation 13 



10 



where i is a divisor of n and /^ J, h^J^ are given in Definition\^ 
For example, 

B^B{5T) = -T2(2e^^ + Se^^ + 3e^ + 2)B(5T) + (-2e^ + 3)Bl 



Relation Uni writes B'^(£T)B(r7,T) as a Q[T, e-^]-linear combination of B(nT) 
and B^+i(£T). Multiply Relation [H by B{nT) and using Relation ^ we 
are able to write B(mT)B^(nT) as a Q[r, e^^] -linear combination of B(nT), 
B^(nT), B^(£T), B^(£T) and B(mT). Repeating the process, we are able to 
write B(r7iT)B'^(nT) as a Q[T, e"^]-linear combination of B(nT), ■ ■ ■ , B^{nT), 
B^(£T), •••, B'^"''-'^(£T) and B(r7iT). Thus we are able to write the product 
B(miT) ■ ■ ■ B(r?7,sT) as a Q[T]-linear combination of elements of the form 
B"(/cT)e""^. Using Relation [TOl the product B(r7iiT) ■ ■ ■B(msT) can be writ- 
ten as a D-linear combination of elements of the form B(6T)e""^. 

Here we carry out the computation for B(2T)B(3T)B(5T). Relations [T2] and 
[13] together with the relation (e-^ — 1)B^ = TB^ (from Relations [7| and [8]) give 
rise to 



B(2T)B(3T)B(5T) 
= B^B{5T) + -T(e^ - l)B(3r)B(5r) - -TB(2r)B(5r) 

= -T^{2e^^ + 3e^^ + 3e^ + 2)B(5T) - 2TB2 + B^ 
5 

+-T{e^ - 1)B2 + -T\e^ - l)(-2e3^ + e^^ - e^ - 3)B(5r) 
3 5 

+ HT2(e2^-2e^ + l)B(3T) 

-^^^B^ - -T\2e'^ - e^^ + e^ - 2)B(5T) + -T^B{2T). 

By Relation [g 



(e^ - l)(-2e^^ + e^^ - e^ - 3)B(5T) = 5(6^^ + 1)B(5T) - lOB, 
(e^^ - 2e^ + l)B(3r) = -3e^B(3T) + 3B. 

Together with the relation [e^ — 1)B^ = TB, we obtain 



B(2T)B(3T)B(5T) = B^ - -TB^ + — T2B(2r) - —T'^e^B{'iT) 

+ V(e3^ + e2^ + 3)B(5T). 





11 



By Relation M, 



B(2T)B(3T)B(5T) = Qt^^ + 5T'-^ " ^^ + l^' - ST + l) B 

4 3 

5 5 5 
In terms of Bernoulli polynomials, 



i,j,k>o V^^^'^J 

i+j+k=n 

]^{n - l){n - 2)Bn + ^n{n - 2)fi„_i 

Related to the representation of the product B(2T)B(3T)B(5T)B(6T), an 
identity of Bernoulli polynomials can be found in [64 Proposition 4] . 



4 Variations 



There are identities of Bernoulli numbers, which do not come directly from 
the algebraic structures of ^. For an example, we replace the underlying field 
Q by the field M of real numbers. 

Relation 14 For s eM., 

B{sT)B{{l-s)T) = (1-.) (b{sT) + ^) B+. [b((1 - s)T) + ^^^) B. 



This relation follows from direct computation and is equivalent to 



(B(sT)-l)(B((l-s)T)-l) 

(1 - .) (B(.r) + ^ - 1) B + s (^B((l - .)r) + ii^ - 1 j B 

+l + B-B(sT)-B((l-s)T). 



12 



Equating the coefficients of T" for n > 0, we obtain 



i+j=n i+2k=n 

+ (1 - s'^ - (1 - sD ^. 

w. 
Dividing tliis identity by s(l — s), tlie integrations 



^ ' s{l-s) (z + j-1)! ' 

give rise to 



.5^0 (--!)! ^ J - ,^^ A;£!(2A:)!+ ^^ nl' 

i+j=n £+2k=n 

For n > 4, multiplying tfie above equation by {n — 1)! and rearranging tlie 
summations, we obtain Miki's identity [IT] 



EBi Bn—i 2 „ 1 .r — , I n\ /„ Bh „ B„—k 
— ^^ = -HnBn + - y , 5„_fc— ^ + Bk—^ 

1=2 ^ ^ ~ ^ ra ra ^ \^Ky V k n — k 

2 „ ■! — > I n\ Bu B„—k 

n ^ \kj k n — k 

Tlie above proof of Miki's identity is taken from M. C Crabb [2J, which is 
essentially a distillation of an argument given by G. V. Dunne and C Schubert 
[1]. The argument works also for the generalization given by I. M. Gessel [^. 

Identities of Bernoulli numbers may come from coefficients in non-trivial ele- 
ments in 03. Consider 



jk jk+i 

^ = {2k + l)T'4- + T'^^''' 



We claim that the coefficient of T^^^ in e-^y^B is always zero. For a proof, we 
use the formula 

to write ip as the sum of ipi and ip2, where ipi = {k + 1)T^ {d^ / dT'') and 
(P2 = T{d/dT)-T^{d^/dT''). Note that the coefficients of T^+^ in {k + l)\e^ipiB 



13 



and {k + l)\e'^ip2B are E [''1j^){k + l)Bk+j and E C'+^^jEfc+j, respectively. 
Thus the identity 

j=o \ 3 / 
discovered by M. Kaneko [9J is equivalent to our claim. 
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